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Abstract. We review the complex differential geometry of the space of ori- 
ented affine lines in M.^ and give a description of Hamilton's characteristic 
functions for reflection in an oriented surface in terms of this geometry. 



1. Introduction 

In a series of classic papers [5] [6] [7] , Hamilton presented his system of geometric 
optics whereby the evolution of systems of rays undergoing reflection and refraction 
at interfaces of homogenous media was described, both in general terms and in 
specific cases. In particular, Hamilton introduced three characteristic functions for 
an optical system by which he could give algebraic expressions for reflection and 
refraction of a given system of rays off simple geometric surfaces, such as spheres 
or other surfaces of revolution. This he carried out in closed form for some systems 
and in approximation for others. For a beautiful exposition of this work, see Synge's 
book [9]. 

The purpose of this paper is to review the complex differential geometry of 
the space of oriented affine lines in R'^ and to give a description of Hamilton's 
characteristic functions for reflection in an oriented surface in terms of this 
geometry. This formalism, developed in [1], has proven useful in the scattering of 
waves off surfaces [2] , as well as the computation of caustics [3] and approximations 
of the Casimir force [4]. 

In the following section we explain the background of this formalism. Section 
3 introduces and describes the three characteristic functions for reflection in an 
arbitrary surface in terms of the complex geometry on the space of lines (Theorems 
3 to 5). 

2. The Complex Geometry of Reflection 

2.1. The Space of Oriented Lines. The fundamental building block of geometric 
optics in a homogeneous isotropic medium is a ray, or oriented straight line, along 
which light is postulated to propagate. The set L of all such oriented lines in 
Euclidean M'^ is a 4-dimensional manifold that can be identified with the total 
space of the tangent bundle to the 2-sphere [8]. To see this, recall that an oriented 
line is uniquely determined by a pair of vectors: its direction vector V and the 
vector U from the origin to the point on the line closest to the origin. The former 
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can be taken to be a unit vector and is perpendicular to the latter. Thus 
L ^ { (C7,f) e X I \V\ = 1 U ■¥ = [>}, 

which is clearly homeomorphic to the tangent bundle to the 2-sphere. 

Our approach to geometric optics is to study optical systems as submanifolds 
of the space L. In order to describe this algebraically we now introduce a local 
coordinate atlas on L. 

Let ^ be the local complex coordinate on the unit 2-sphere in M.'^ obtained by 
stereographic projection from the south pole. In terms of the standard spherical 
polar angles {9, (j)) about the north pole, we have ^ = tan(|)e*'^. 

This can be extended to complex coordinates rj) on L minus the tangent space 
over the south pole, as follows. Note that a tangent vector X to the 2-sphere can 
always be expressed as a linear combination: 

^ d a 

for a complex number rj. We can thus identify the real tangent vector X on the 
2-sphere (and hence the ray in M.^) with the two complex numbers (^,r/). Loosely 
speaking, ^ determines the direction of the ray, and ij determines its perpendicular 
distance vector to the origin. 

The coordinates r]) do not cover all of L - they omit all of the lines pointing 
directly downwards. However, the construction can also be carried out using stere- 
ographic projection from the north pole, yielding a coordinate system that covers 
all of L except for the lines pointing directly upwards. The whole of the space of 
oriented lines is covered by these two coordinate patches and the transition func- 
tions between them are holomorphic. In what follows we work in the patch that 
omits the south direction. 

Definition 1. The canonical bundle map tt : L ^ is given by mapping an 
oriented line to its direction. Here, and throughout, we write for with the 
standard complex structure, and denote the composition of projections L x R ^ 
L Pi by TTi. 

One of the key relationships between L and is the map $ : L x M — > M^. 

Definition 2. The map $ takes ((^, 77), r) G L x R to the point on the oriented line 
{£,, rj) in that lies an affine parameter distance r from the point on , rj) closest 
to the origin. 

If $((^,7/),r) = {z{^,ri,r),t{^,r],r)), where z = x'^ +ix'^, t = and {x^, x"^, x^) 
are Euclidean coordinates in R^, then we have the following coordinate expressions 
[1]: 



(2.1) 



This map is of crucial importance when describing surfaces in R^, as we explain 
below. It is also useful to invert these relationships: 
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These give the perpendicular distance vector from the origin and the distance r for 
a hnc through the point (z, t) E M'^ with direction ^ G P"'^. 
For later use we introduce the following definition: 

Definition 3. For a subset iiT C M^, the subset C L is defined to be the set of 
oriented lines that intersect K: 

Ua' = { 7 e L I $(7, ■) n is: 7^ } . 

Example 1. For p G R'^ the set Up is the set of oriented lines through p, which is 
a sphere in L. In terms of the coordinates above, if p is then 

U(,,t) = { (^,77) G L h = i (z - 2tC - }■ 

Example 2. For a surface S" C M"^, Us is a 4-dimensional subset of L whose 
boundary is 

dVs = { 7 G L I $(7, •) tangent to S }. 

In its simplest form, an optical system is a 2-parameter family of rays. We refer 
to a 2-parameter family of oriented lines in M'^ as a line congruence, which, from 
our perspective, is a surface E in L. For example, a point source corresponds to 
the 2-parameter family of oriented lines that contain the source point, which thus 
defines a 2-sphere in L - the sphere in Example 1. 

For computational purposes, we now give explicit local parameterizations of the 
line congruence. In practice, this will be given locally by a map C — > L : ^ 
(^(/i, /2), 77(^, /i)). A convenient choice of parameterisation will often depend upon 
the specifics of the situation, but our formalism holds for arbitrary parameterisa- 
tions. 

We now describe how to construct surfaces in using line congruences. Given 
a line congruence E C L, a map r : E ^ M determines a map E ^ by (^, rj) 
$((^, 77), r(^, 77)) for (^,77) G E. With a local parameterization n of E, composition 
with the above map yields a map C — > which comes from substituting ^ = 
^{fi,fl), r/ = ri{p,fl) and r — r{fi,fi) in equations (2.1). In other words, we choose 
one point on each line of the congruence, thus forming a surface. 
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Definition 4. An oriented surface S in E'^ gives rise to the map : 5* ^ LxK 
which takes a point p G S" to the oriented normal hne to S crossed with the distance 
of p to the point on the normal lying closest to the origin. 

The dual picture of light propagation is to consider the wavefronts, or surfaces 
that are orthogonal to a given set of rays. However, not every line congruence has 
such orthogonal surfaces - indeed, most don't. We refer to line congruences for 
which orthogonal surfaces exist as normal congruences. 

Theorem 1. [1] A parameterised line congruence (^(fi, p,),ri{ii, ft)) is orthogonal to 
a surface in iff there exists a real function r{fi, p.) satisfying: 

where d — 

If there exists one solution, there exists a 1 -parameter family generated by the 
real constant of integration. An explicit parameterization of these surfaces in M.^ is 
given by inserting {£^{^, fl),ri{^, p,)) and r = r{fj,, p,) in (2.1). 

Given a surface S with coordinates ^ = S^q{^,P), rj — 7]o{fi,p) and r = ro{iJ,,p), 
the set Vs is [2]: 

(2.4) 

2.2. Reflection. Given an oriented surface S* in R"^, reflection in S determines 
a map TZs : Us Ug. 




The algebraic expression for this map is: 

Theorem 2. [2] Let S be an oriented surface S in R'^ with parameterised normal 
line congruence ^ = ^o(MjM); V — Voilijp-) o,nd r = ro(/i,/l) satisfying (2.3). Then 



ON HAMILTON'S CHARACTERISTIC FUNCTIONS FOR REFLECTION 



5 



TZs -Vs Us takes (^i,rii) to TZs{S,i,rii) = (67?72)- 

(2.6) 

where {£,o,T]o,rQ) e ^s^{S). 

If ri and r2 are the distances of the point of reflection from the point closest to 
the origin on the incoming and reflected rays, respectively, then 

3. Hamilton's Characteristic Functions For Reflection in a Surface 

Consider a ray (Cii'7i) entering an optical instrument and emerging, after a 
number of reflections and refractions, as the ray (^2, '72)- Let pi and p2 be arbitrary 
points on the incoming and outgoing rays, qi and 52 be the points on the incoming 
and outgoing rays that he closest to the origin (respectively). 




Hamilton defined three different functions, which he termed characteristic func- 
tions and denoted T, V and W . 

Definition 5. The function T, which is referred to as the angle characteristic 
function, takes the directions ^1 and ^2 to the distance along the rays between qi 
and q2. 

Definition 6. The function W , which is referred to as the mixed characteristic 
function, takes the point pi and the direction ^2 to the distance along the rays 
between pi and (72 ■ 

Definition 7. The function V , which is referred to as the point characteristic 
function, takes the points pi and p2 to the distance along the rays between pi and 

P2- 
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The function V is, of course, commonly called the Hamiltonian of the optical 
system. 

In the above, while the points are on the initial or final rays, they may be on 
the opposite side of the initial or final interface and so lie at points which the light 
may never physically pass through. However, as long as the associated distances 
are counted negatively, the functions are well defined and have all of the formal 
properties required for our purposes. 

The exact domain of definition of these functions will depend upon the optical 
instrument under consideration. In addition, the characteristic functions may be 
multivalued, since there may be more than one sequence of rays with the properties 
required. 

The domain of T is obviously a subset of x , while the domain of V" is a 
subset of X P-^ and the domain of is a subset of x K^. More specifically, 
for reflection in a surface: 

Proposition 1. Let S be an oriented immersed surface in K^, then for reflection 
in S: 

DomT={ (^1,6) e X Pi \TZs (7r-i(Ci) n Us) n Tf-ife) ^ } 
Dom W= { eM' xPi I 7es(Up, nUs)n7r-i(C2) ^0 } 

Dom V= { {pi,P2) G x r3 I 7^_5 (jj^^ n Us) n Up, ^ } . 

Proof. A pair of directions ^1,^2 £ P^ li<^ in the domain of T iff there exists a line 
with direction ^1 which, when reflected in 5", has direction ^2- Now 7r~"'^(^i) is the 
set of oriented lines with direction ^1, while 7r^^(^i) flUs is the set of oriented lines 
with direction ^1 that intersect 5*. Thus, reflecting this set gives TZs H Us) 

and, for (^1, ^2) to be in the domain of T, the reflected set must contain a line with 
direction ^2- In other words, TZs ("""^^(^i) Us) n 7r~^(^2) 7^ 0, as claimed. 

For W, {pi , ^2) G X P^ is in the domain of W iff there exists a line through pi 
whose reflection in S has direction ^2 • The oriented lines through pi that intersect 
S are Up^ n Us. Thus, for (pi, ^2) to be in the domain of W, the reflection of this 
set, TZs (Upi nUs) must contain a line with direction ^2- That is, we must have 
TZs (Upi n Us) n 7r-i(^2) ^ 0, as claimed. 

Finally, the set of oriented lines through pi that intersect S is Up^ D Us- Thus, 
for (pi,P2) to be in the domain of V, the reflection of this set must contain an 
oriented line that passes through p2- That is, TZs (Up^ fl Us) n Up, 7^ 0. □ 

The following propositions describe the domains of the characteristic functions 
in terms of our coordinates. We assume that the oriented surface S has normal 
congruence parameterised by /i ^ (Co(Mj P), '7o(a'-, ft)) and r = ro{n, j2). 

Proposition 2. Given $1,^2 G P\ define ^0 G P^ by 

. ^ - 66 + la - 61 [(1 + ac"i)(i + 66)]^ 
a(i+66)-6(i+a6) ■ 

Then the domain ofT is given by: 

Dom T= { (a, 6) e Pi X Pi I ^0 e TTi o$-i(5) } . (3.2) 

Proof. Given two directions 6:6 G P"'^ , the direction ^ through which ^ must 
be reflected to obtain ^ (or vice versa) is found by solving equation (2.5) for ^o- 
The result is (3.1), and the question then reduces to whether the normal to S 
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has direction at any point. Now, the map tti o : 5 ^ takes a point 
on an oriented C"^ surface S to the direction of the normal at the point (i.e. the 
classical Gauss map of S). Thus, for (^i, ^2) to be in the domain of T we must have 

eoe7rio$-i(^). □ 

Proposition 3. The domain ofW is obtained as follows. Given pi = {zi,ti) G 
and ^2 G P^, solve (if possible) the single complex equation for /j.: 

zi[{l - ^o^o)C2 - 2^0]' - 2ti[(l - eo^o)^2 - 2eo][2^o^2 + 1 - ^o^o] 

-^i[2^o6 + l-Coeo]' 

= 2(1 + ^oCo) ((fo - 6)'% -(1 + ^o6)'% + (fo - 6)(1 + eo6)(l + ^ofo)ro) • 

(3.3) 

Then {pi,^2) is in the domain ofW iff^o{fJ-,p-) G tti o $^^(5). 

Proof. An oriented line (^i,r/i) passes through the point {zi,ti) iff {of. equation 
(2.2)) 

= i (zi - 2tiei - zi^?) , 
while inverting the reflection law (2.5) gives 

^ 2^0^2 + 1-^0^0 

(l-W6-2Co" 
For the incoming ray to intersect the surface S we must have {of. (2.4)) 

Thus, for an oriented line through pi to have reflected direction ^2, the equation to 
be solved is obtained by substituting the first two equations above in the last one. 
The result is as claimed. □ 

Proposition 4. The domain of V is determined as follows. Given two points 
Pi = {zi,ti) and p2 = {z2,t2) G 'S? , solve (if possible) the following two complex 
equations for n and ^i: 

Z2[{1 - CoCo)a - 2Co]' - 2t2[(l - Co^"o)6 - 2Co][2^o6 + 1 - ^o^o] 
+^2[2Coa + 1 - ^oCo]' = -(1 + ^o^o)'(^i - 2*1^1 - ziC?) 

+ 4(fo - a)(l + ^o^"i)(l + €o€"o)r-o (3.4) 

and 

(1 + €0^0)' {zi - 2ii6 + ^1^1') = 2(1 + ^o^ifvo - m - Cifm 

-2(l + ^oa)(€o-6)(l+CoC"o)ro. 

Then (pi,P2) is in the domain of V iff ^q{ij,, fi) £ tti o^g^[S). 
Proof. As in the previous proposition, we must have 

V'i — 2 (^^ 2i^^2 Zi^^ ) , 

this time for i = 1,2. Substituting these into equations (2.4) and (2.6), with the 
aid of (2.5), yields the two equations in the proposition. □ 



(3.5) 
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We now compute the characteristic functions for reflection in an oriented 
surface S. Assume that the normal line congruence to S is parameterized by ^ = 
^o(/«, P), V = A) and r = ro(/x, p). 

The angle characteristic function for reflection in a surface S is given by: 

Theorem 3. Hamilton's angle characteristic function for reflection in S is 
where is given by (3.1) and (^o,'?o,?'o) S ^s^('^)- 

Proof. Given and ^2, from the reflection law (2.5) we solve for and get the 
direction of the normal to the surface at the point of reflection - equation (3.1). 
It is possible that the surface S never has this normal direction, in which case the 
function T is not deflncd at (^1,^2)- On the other hand, there may be more than one 
point on the surface with normal direction in which case T will be multivalued. 
Assume that it is defined and that the point of reflection is {zo,to) which lies at 
r = ro along the normal ray (^cf^o)- 

As before, let ri and 7-2 be the distances of the point of reflection from the 
point closest to the origin on the incoming and reflected rays, respectively. Then 
T(a,6)=±|ri-r2|andby(2.7) 

2(|i+c'o6p-igo-6n ^ 
(i+Mo)(i+eiei) 

Substituting equation (3.1) in this yields the result after some simplification. □ 



On the other hand the mixed characteristic function for reflection is determined 
by: 

Theorem 4. Hamilton's mixed characteristic function for reflection in S is: 

W^((2i,ii),6) = 



eizi + eizi + (i-aa)ti ^ 2ia-6i 

H 1 1 — r^'o 



1 + 66 [(i + ^^^-)(i + ^26) 



where 



^^=(1-^0^0)6-26' ^'-'^ 
and ^0 G "^1 ° ^^^{S) is a solution of (3.3). 

Proof. From Proposition 3, (2:1,^1), ^2) is in the domain of W iff there exists a 
direction ^0 G tti o $^^(5') that solves equation (3.3). Assuming that such exists, 
equation (3.6) follows from inverting the reflection law (2.5). 
It is clear that 

W^((2i,ii),6) = ±|si + n -7-2!, 

where si is the distance from pi to q\. Prom the second equation of (2.2), we have 
that 

6-21 + 6^1 + (1-66)^1 

= rT66 ' 
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while from the proof of the previous theorem, we have that 

2ia-6i 

ri-r2 = z : — jTo- 

The result follows. □ 



Finally, the point characteristic function for reflection can be computed using; 
Theorem 5. Hamilton's point characteristic function for reflection in S is: 



where 



2|ei-6i 

si - S2 H z 1 — r'^o 



^iZi + ^iZi + (1 - ^i^i)ti 



(3.7) 



1 + 

fori = 1,2, 

, 2^0^ + 1 -^0^0 
(l-^o^o)6-2^"o' 

and ^0 and solve equations (3.4) and (3.5), and S tti o $^^(5). 

Proof. From Proposition 4, (zi, ti), (z2, ^2)) is in the domain of V iff there exist 
directions ^0 and ^1 that solve equations (3.4) and (3.5), and ^0 G tJ"! ° 'J's^('S'). 
Assuming that such exists, equation (3.8 is just the reflection law (2.5), and 

V{{ziM),{z2,t2)) = - S2 + ri - r2\ . 
Finally, the distance from pi to Qi is given by (3.7). □ 
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